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Abstract. We establish an isoperimetric inequality for the fundamental tone 
(first nonzero eigenvalue) of the free plate of a given area, proving the ball 
is maximal. Given r > 0, the free plate eigenvalues u> and eigenfunctions 
u are determined by the equation AAu — tAu = uju together with certain 
natural boundary conditions. The boundary conditions are complicated but 
arise naturally from the plate Rayleigh quotient, which contains a Hessian 
squared term |D 2 u| 2 . 

We adapt Weinberger's method from the corresponding free membrane 
problem, taking the fundamental modes of the unit ball as trial functions. 
These solutions are a linear combination of Bessel and modified Bessel func- 
tions. 



1. Introduction 

Laplacian and bi-Laplace operators are used to model many physical situations, 
with their eigenvalues representing quantities such as energy and frequency. The 
eigenvalues fj, of the Neumann Laplacian on a region fl determine the frequencies 
of vibration of a free membrane with that shape. If f2* is the ball of same volume 
as f2, then we have 

< Mi(^*) with equality if and only if Q is a ball. 

First conjectured by Kornhauser and Stakgold [13], this isoperimetric inequality 
was proved for simply connected domains in M 2 by Szego 26, 25] and extended to 
all domains and dimensions by Weinberger |32) . 

The main goal of this paper is to establish the analogous result for the eigenvalues 
u! the free plate under tension. That is, of all regions Q with the same volume, we 
have the bound 

Wi(f2) < ll>i(Q*) with equality if and only if Q is a ball. 

Our proof relies on the variational characterization of eigenvalues with suitable trial 
functions. Similar to Weinberger's approach for the free membrane, we take our 
trial functions to be extensions of the fundamental mode of the unit ball. However, 
because the plate equation is fourth order, finding the trial functions and estab- 
lishing the appropriate monotonicities is significantly more complicated than in the 
membrane case. The eigenmodes of the unit ball are identified in our companion 
paper [5J, where we also establish several properties of these functions that are 
used in the proof of the isoperimetric inequality. If the reader is satisfied with a 
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numerical demonstration, the needed properties of ultrasphcrical Bessel functions 
can be verified in any given dimension using Mathematica or Maple. 

The boundary conditions of the free plate are not imposed, but instead arise 
naturally from the Rayleigh quotient. It is therefore extremely important that we 
begin with the correct Rayleigh quotient for the plate, which includes a Hessian 
term. These natural boundary conditions have long been known in the case d = 2 
(see, eg, |33|); we include in this paper their derivation for all dimensions. 

The fundamental tone of the ball is extremal for other physically meaningful 
plate boundary conditions. The ball provides a lower bound for the clamped plate 
eigenvalues [23 HH1 EH E] • The methods used by Talenti, Nadirashvilli, Ashbaugh 
and Benguria to prove the clamped plate isoperimetric inequality are quite differ- 
ent than those for the free plate and membrane and only establish the bound in 
dimensions 2 and 3. The problem remains open for dimensions four and higher, 
with a partial result by Ashbaugh and Laugesen [3] . For an overview of work on 
the clamped plate problem, see 10, Chapter 11, p. 169-174] and [HI p. 105-116]. 

Other plate boundary conditions include the simply supported plate, hinged 
plate, and Neumann boundary conditions. Plate problems are fourth-order and 
generally more difficult than their second-order membrane counterparts, because 
the theory of the bi-Laplace operator is not as well understood as the theory of the 
Laplacian. Verchota recently established the solvability of the biharmonic Neumann 
problem |31j ; these boundary conditions arise from the zero-tension plate and allow 
consideration of Poisson's ratio, a measure of a material property that we take to 
be zero for our free plate. Supported plate work includes Payne [23 and Licar and 
Warner |15j . who examine domain dependence of plate eigenvalues. It would be 
natural to conjecture an isoperimetric inequality for the simply supported plate, 
although there does not seem to be any work on this problem to the best of our 
knowledge. Work with hinged plates includes Nazarov and Sweers [20] . 

Other notable mathematical work on plates includes Kawohl, Levine, and Velte 
[TT] , who investigated the sums of low eigenvalues for the clamped plate under ten- 
sion and compression, and Payne |23| , who considered both vibrating and buckling 
free and clamped plates and established inequalities bounding plate eigenvalues by 
their (free or fixed) membrane counterparts. For a broad survey of results, see 

mm- 



2. Formulating the problem 

We now develop the mathematical formulation of the free plate isoperimetric 
problem. Let fl be a smoothly bounded region in M. d , d > 2, and fix a parameter 
t > 0. The "plate" Rayleigh quotient is 

S n \D*u\ 2 +r\Du\*dx 

« QM = ' 

Here \D 2 u\ = {J2jk u ljx k ) 1/2 is the Hilbert- Schmidt norm of the Hessian matrix 
D 2 u of u, and Du denotes the gradient vector. 

Physically, when d = 2 the region f2 is the shape of a homogeneous, isotropic 
plate. The parameter r represents the ratio of lateral tension to flexural rigidity of 
the plate; for brevity we refer to r as the tension parameter. Positive r corresponds 
to a plate under tension, while taking r negative would give us a plate under 
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compression. The function u describes a transverse vibrational mode of the plate, 
and the Rayleigh quotient Q[u] gives the bending energy of the plate. 

From the Rayleigh quotient 0, we will derive the partial differential equation 
and boundary conditions governing the vibrational modes of a free plate. The 
critical points of ([!]) are the eigenstates for the plate satisfying the free boundary 
conditions and the critical values are the corresponding eigenvalues. The equation 
is: 

(2) AAm - tAu = lou, 

where ui is the eigenvalue, with the natural {i.e., unconstrained or "free") boundary 
conditions on dtt: 

d 2 u 

(3) Mu := -—r = 

on 2 

(4) Vu := r 9 £ - dW an (Pan [(D*u)n\) - ^ = 

Here n is the outward unit normal to the boundary and div^n and gradg^ are the 
surface divergence and gradient. The operator Pqq projects onto the space tangent 
to 90. 

We will prove in a later section that the spectrum of the Rayleigh quotient Q is 
discrete, consisting entirely of eigenvalues with finite multiplicity: 

= LOq < Ul\ < Ll>2 <•••—> oo. 

We also have a complete L 2 -orthonormal set of eigenfunctions uq = const, Ui, 112, 
and so forth. 

We call Mi the fundamental mode and the eigenvalue lj\ the fundamental tone; 
the latter can be expressed using the Rayleigh-Ritz variational formula: 

Wi(0) = min{Q[u] : u G H 2 (ti), / udx = 0}. 

In general, the kth eigenvalue is the minimum of Q[u] over the space of all functions 
u i 2 -orthogonal to the eigenfunctions uq, u\ r . . , u^-i- Because uq is the constant 
function, the condition u _L uq can be written J^udx — 0. Note that in the limiting 
case t = 0, the first d+ 1 eigenvalues of O are trivial because Q[u] =0 for all linear 
functions u. We therefore need the tension parameter r to be positive in order to 
have a nontrivial inequality. 

The eigenvalue equation ^ can also be obtained by separating the plate wave 
equation 

(j> tt = —AAcf) + t A</>, 

by the separation <fi(x, t) — u(x) cos(^/uJt). The eigenvalue to is therefore the square 
of the frequency of vibration of the plate. 

3. Main result 

The main goal of this paper is to prove an isoperimetric inequality for the funda- 
mental tone of the free plate under tension. Let O* denote the ball with the same 
volume as our region f2. 

Theorem 1. For all smoothly bounded regions of a fixed volume, the fundamental 
tone of the free plate with a given positive tension is maximal for a ball. That is, 
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Figure 1. The fundamental tone wi(fT) of the disk (middle 
curve) and of another region Wi(f2) (bottom curve). The dashed 
straight line is tangent to wi(J7*) at r = 0. 



if t > then the first nonzero eigenvalue lj of AAu — tAu — luu subject to the 
natural boundary conditions ^ and Q satisfies 

(5) lui(Q) < wi(0*), with equality if and only if is a ball. 

The proof of Theorem [T] will proceed from a series of lemmas, following roughly 
this outline. A more detailed summary follows. 

• Section 8 - Defining the trial functions, showing concavity of the radial 
part of the trial function, and evaluating the Rayleigh quotient 

• Section 9 - Proving partial monotonicity of the Rayleigh quotient 

• Section 10 - Establishing rescaling and rearrangement results, and proving 
the theorem. 



Adapting Weinberger's approach for the membrane [32], we construct in Lemma 13 
trial functions with radial part p matching the radial part of the fundamental mode 



of the ball. We follow by proving in Lemma 14 a concavity property of p that 



will be needed later on. We next bound the eigenvalue w by a quotient of inte- 



grals over our region f2, both of whose integrands are radial functions (Lemma 161 



These integrands will be shown to have a "partial monotonicity". The denomina- 



tor's integrand is increasing by Lemma 17 and the numerator's integrand satisfies 



a decreasing partial monotonicity condition by Lemma 18 

The proof of Lemma [18] becomes rather involved and so is contained in its own 
section and broken into two cases, Lemma [T9] for large r values, and Lemma [20] for 
small values of r. The latter in turn requires some facts about particular polynomi- 



als, proved in Lemmas 21 and 22 We then exploit partial monotonicity to see that 
the quotient of integrals is bounded above by the quotient of the same integrals 
taken over fT, by Lemma [23] Finally, we conclude that the quotient of integrals 
on fl* is in fact equal to the eigenvalue lj* of the unit ball. From there we deduce 
the theorem. 
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4. Existence of the spectrum and regularity of solutions 

Our first task is to investigate the spectrum of the fourth-order operator asso- 
ciated with our Rayleigh quotient Q in ([I]). In this section we show the spectrum 
is entirely discrete, with an associated weak eigenbasis. We will then establish reg- 
ularity of the eigenfunctions up to the boundary and derive the natural boundary 
conditions. 

For this section only we will allow r to be any real number. We continue to 
require C R d to be smoothly bounded unless otherwise stated. 

The existence of the spectrum. We consider the sesquilinear form 

r. d 

a(u,v) = / u XiXi v XiX . + t(Du ■ Dv) dx 

in L 2 (f2) with form domain H 2 (£l). Note the plate Rayleigh quotient Q can be 
written in terms of a, with Q[u] = a(u, u)/||it|| 2 2 ■ 

Proposition 2. The spectrum of the operator A 2 — rA associated with the form 
a(-, •) consists entirely of isolated eigenvalues of finite multiplicity uj < u)\ < w 2 < 
• • • — > oo. There exists an associated set of real-valued weak eigenfunctions which 
is an orthonormal basis for L 2 (Q). 

Proof. By Cauchy-Schwarz, the form a(-, •) is bounded on H 2 (Q), and so is contin- 
uous. We will show the quadratic form a(u, u) is coercive; that is, for some positive 
constants c\ and c 2 a(u,u) + ci||u|| 2 > C2||m|| 2 ^2(q). By the boundedness of a on 
H 2 , this is equivalent to showing the norm associated with a, 

IMIa = a(u,u) +ci||u|| 2 
is equivalent to || • |||p(Q), and hence a is closed on H 2 (Q). Because f2 is smoothly 
bounded, H^Q) and L 2 (fi) can be extended to H 2 (R d ) and L^R d ) respectively. The 
space H 2 (R d ) is compactly embedded in L^R d ). Then by a standard result (see 
e.g., Corollary 7.D [23 p. 78]), the form a has a set of weak eigenfunctions which 
is an orthonormal basis for L 2 (J7), and the corresponding eigenvalues are of finite 
multiplicity and satisfy 

(6) u>i < ll>2 < • • • < w tt — > 00 as n — > 00. 
For r > 0, coercivity of the form a is easily proved: 

a(u,u) +t\\u\\ 2 > \\D 2 u\\ 2 + t\\Du\\ 2 + t\\u\\ 2 
> min(r, l)||«||^a, 

where all unlabeled norms are L? norms on fi. 

To prove coercivity when r < 0, we must somehow arrive at a positive constant 
in front of the \Du\ 2 term. We cannot use Poincare's inequality on the \D 2 u\ term 
as this will introduce terms involving the average value of Du. Instead, we will 
exploit an interpolation inequality. 

By Theorem 7.28 of [HI p. 173], we have that for any index 1 < j < n and any 
£ > 0, 

(7) \\d Xj u\\ 2 LHn) < e\\uf H2(n) + Ce^WuWl^ 
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with C = C(f2) a constant. Replacing e by e/d and summing over j, we see 
\\D 2 u\\h > [I - l) \\Du\\% 2 -(^ + l 

Fix 5 e (0, 1). Let K > 0. Then 

a(u,u)+K\\u\\ 2 L2 = \\D 2 u\\h - \t\\\Du\\ 2 l2 + K\\u\\l 2 

S , , A „„ „ 2 , / ^ Cr) , > 



> (1 - <5)||^||i 2 + I - - 5 - |r| j ||Dti||* a + Iff - _ - * l|| u 

We can choose our e small and our K large so that the minimum is positive, which 
proves coercivity. For example, for 5 = 1/2, we need to take e < 1/(1 + 2|r|) and 

K> §(c+1 + 2|t|). Thus a is coercive for all r. 

Now suppose u is a weak eigenfunction corresponding to eigenvalue uj. Because 
u is real-valued, by taking the complex conjugate of the weak eigenvalue equation 
we see that u is also a weak eigenfunction with the same eigenvalue. Thus the real 
and imaginary parts of u are both eigenf unctions associated with w, and we may 
choose our eigenfunctions to be real-valued. □ 

Note that for any bounded region f2 and all real values of t, the constant function 
solves the weak eigenvalue equation with eigenvalue zero. For all nonnegative values 
of t, the Rayleigh quotient is nonnegative for all functions and so = wq < w\ < 

■ ■ ■ < Wk < When r = 0, the coordinate functions xi, . . . ,Xd are also solutions 

with eigenvalue zero, and so the lowest eigenvalue is at least d + 1-fold degenerate, 
as noted in the introduction. Taking instead r > 0, the Raleigh quotient shows 
that the fundamental tone uji is positive, and so we have: 

= ujq < u>i < uj2 < • • • < U) n — » oo as n — > oo. 

Regularity. We aim to establish regularity of the weak eigenfunctions by appeal- 
ing to interior and boundary regularity theory for elliptic operators. 

Proposition 3. For any r£l and smoothly bounded fl, the weak eigenfunctions 
of A A — rA are smooth on fi. 

Proof. Let u be a weak eigenfunction of A with associated eigenvalue w; by Propo- 
sition^ we have u € T(a) = H 2 (fl). Then by a theorem in J2TJ p 668], we have 
u G ZP(0) for every positive integer k. Thus we have u £ H k (Q) for all k e Z + , 
and so u £ C°°(n). 

Regularity on the boundary follows from global interior regularity and the Trace 
Theorem (see, for example, [301 Prop 4.3, p. 286 and Prop 4.5, p. 287.]). Thus we 
have u £ C°°(f2), as desired. □ 

5. The Natural Boundary Conditions 

In this section, our goal is to derive the form of the natural boundary conditions 
necessarily satisfied by all eigenfunctions. Consider the weak eigenvalue equation 
for eigenfunction u with eigenvalue oj and some test function c/> g C^°(fl): 

n ij 
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Because the eigcnfunction u is smooth, we may use integration by parts to move 
most of the derivatives on (j) to u; this gives us a volume integral and two surface 
integrals that must vanish for all <fi. 

The natural boundary conditions are rather complicated in higher dimensions, 
and so we state the two-dimensional case first. The boundary conditions in this 
case have been known for some time: see, for example, |33j 

Proposition 4. (Two dimensions) For C M 2 , the natural boundary conditions 
for eigenf unctions of the free plate under tension have the form 

d 2 u 

Mu := — -= = 

Vu . = _ d ( Au ) _ _d_ ( Q2u _ K{s) dv ^\ = 
dn dn ds \dsdn ds J 

where n denotes the outward unit normal derivative, s the arclength, and K the 
curvature of dQ . 

We also look at one example of the natural boundary conditions for a region 
with corners. Notice that an additional condition arises at the corners! 

Proposition 5. (Rectangular region in two dimensions) When £1 C M 2 is a rect- 
angular region with edges parallel to the coordinate axes, the natural boundary con- 
ditions for eigenf unctions of the free plate under tension have the form 

d 2 u , , 

— — ; - =0 at each edge 
dn z 

du d 3 u d(Au) 

T a 7T?~H s =0 on each edge 

on us on on 

Uxy =0 at each corner 

where n and s indicate the normal and tangent directions. 

Finally, we state the natural boundary conditions for a smoothly-bounded region 
in higher dimensions: 

Proposition 6. (General) For any smoothly bounded Q, the natural boundary con- 
ditions for eigenf unctions of the free plate under tension have the form 

Mu := — % = on dn, 

on z 

Vu := t tt~ - divan (Pon [{D 2 u)n} ) - ^ = on dn, 

On V L 1 J On 

where n denotes the normal derivative and divgn is the surface divergence. The 
projection Pqq projects a vector v at a point x on dfl into the tangent space of dn 
at x. 

Proof of Proposition^ Our eigenfunctions u are smooth on n by Proposition 2 and 
satisfy the weak eigenvalue equation a(u, (f>) — ui(u, (j))L 2 (Q) — for all <f> G H 2 (n). 
That is, 



/ 



^ u XiXj <f> XiXj + tD4> ■ Du - cou(j> | dx = 0. 
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As in the membrane case, we make much use of integration by parts. Let n 
denote the outward unit normal to the surface dCl. To simplify our calculations, 
we consider each term separately. 

The gradient term only needs one use of integration by parts: 

/ Du-D4>dx= [ 4>^dS- [ 4>{Au)dx. 
Jn Jan °n J n 

The Hessian term becomes: 



D<\> ■ ((D 2 u)n) - <P^^j dS + j (A 2 u) 



idx. 



I an 

after integrating by parts twice. 

We wish to transform the term involving Dip in the above surface integral us- 
ing integration by parts. Because we are on d£l, we must treat the normal and 
tangential components separately. We can then use the Divergence theorem for 
integration on dQ. 

We note that the surface gradient grad^n equals D — nd n when applied to a 
function (like <p) that is defined on a neighborhood of the boundary. Thus grad asl 
gives the tangential part of the Euclidean gradient vector. Hence, 



an 



D(f>- {[D 2 u)n^j dS 

= J m (i^+^an<p) ■ (n^+Pan [{D 2 u)n}^ dS 

= I 7T ft + (s rad a>> Pan [(D 2 u)n] ) dS 
Jan " n u n \ L 1 i an 

= I 7T ft - 0div a n (P an [(D 2 u)n\) dS, 

by the Divergence Theorem on the surface 9f2. Here (•, -}gn denotes the inner 
product on the tangent space to <9f2. Recall Pqh projects a vector at a point x on 
dQ, onto the tangent space of dft at x. 

Thus for u an eigenfunction associated with eigenvalue lo, we see 

0=1 4>[ A 2 u — tAu — lou ) dx 



d 2 u , / du dAu 



an 



t- 



dn dn 2 \ dn dn 



div 9n (Pan [{D 2 u)n\ ) j dS 



As in the membrane case, this identity must hold for all <p E H 2 (il). If we take 
any compactly supported 0, then the volume integral must vanish; because </> is 
arbitrary, we must therefore have A 2 u — tAu — lou = everywhere. Similarly, the 
terms multiplied by <f> and d(f>/dn must vanish on the boundary. Collecting these 
results, we obtain the eigenvalue equation ^ and natural boundary conditions of 
Proposition [6j □ 

Proof of Proposition^ Here d = 2; take rectangular coordinates (x,y). We param- 
etrize Oil by arclength s and define coordinates (n, s), with n the normal distance 
from dtt, taken to be positive outside f2. Write n(s) and i(s) for the outward unit 
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normal and unit tangent vectors to the boundary. Then Pqq [/in + f$b\ — fit and 
the operators div^n an d grad af2 both simply take the derivative with respect to 
arclength s. That is, for a scalar function f(s), and taking t(s) to be the tangent 
vector to the surface, we have 



grad an /(s) = f'(s) and 
and so we may write 

divan [Pan [(D 2 u)n 



divan(/(s)*(s)) = /'(a). 



a_ 

ds 



t T (D 2 u)n. 



The tangent line to dfl at the point (0, s) in our new coordinates forms an angle 
a = a(s) with the x-axis (see [33J P- 230]); the curvature of dtt is given by K(s) — 
a'(s). Then in rectangular coordinates, the unit tangent vector is (cos a, sin a), and 
the outward unit normal is (sin a, — cos a). Thus we have 

—t T (D 2 u)s = d s ^ sin a cos a(u xx — u yy ) + (sin 2 a — cos 2 a)u xy ^ . 
p. 233], on dQ under our change of coordinates, we have 



By [3 

V'XX 
Uyy 

u 



sm a 



u ss cos 2 a + 2u ns sin a cos a + Ku n cos 2 a — 2Ku s sin a cos a 



u n n cos 2 a + u ss sin 2 a — 2u ns sin a cos a + Ku n sin 2 a + 2Ku s sin a cos a 



L xy 



u nn cos a sin a + u ss cos a sin a + u ns (sin 2 a — cos 2 a) 
+ i^u n cos a sin a — i^Us (sin 2 a — cos 2 a) . 
So after simplification, 



ason as 



This together with the results of Proposition [6] yields the form of Vu given in 
Proposition |4j Mu is unchanged, and so this completes the proof. □ 

Proof of Proposition^ Our previous findings do not completely apply because dfl 
has corners, although our argument proceeds similarly. For convenience of notation, 
we will take fi to be the square [0, l] 2 . 

The Hessian term gives us a condition at the corners. In particular, after inte- 
grating by parts twice, we have: 



T^xx^xx 2li xy 4^xy ^yy^yy dA 



fi^^xxxx ~\~ 2Uxxyy ~l~ ^yyyy^ dA 
{^xx^x 



^xxxfi ~t~ U x y(f)y U xyy (f) 



+ / I Uyy(j) V - Uyyy(j) + U X y<j) ~ 



x=l 
x=0 

0=1 

y=o 



dy 
dx. 



Since 



u xy 4> y dy = u xy <p 



0=1 



0=0 Jo 



Uxyy^dy 
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and 



we obtain 



Jo 



^xy^x dx — U X y(j 



x=l r 1 
I 

*=0 Jo 



1 I 7 1} 1 



idx 



Hxx^xx ~t~ 



+ U V 



,dA 



I A' 



+ 2u 



xxyy 



+ U 



dA 



+ 



{^xx^x fii^^xyy ~t~ ^zxa;)^ 
J i^ u yy ( t ) V ~ 4>{^- u xxy + Uyyy)^ 

y=i 



x=l 
x=0 

y=i 



dy 

dx 



~t~ ^Uxy 



x=l 



x=Q 



y=o 



Because the Divergence Theorem does apply to regions with piecewise-smooth 
boundaries, the gradient term is the same as in the smooth-boundary case. The 
final term above is the only term that depends only on the behavior of u and <p 
at the corners; arguing as before, we obtain the eigenvalue equation and natural 
boundary conditions, with the additional condition 



= U„ 



l 

x=0 



1 

y=o 



That is, we must have u xy = at the corners. 



□ 



Example: natural boundary conditions on the ball 
When SI is a ball, we can simplify the general boundary conditions. 

Proposition 7. (Ball) The natural boundary conditions in the case ft = M d (R), 
the ball of radius R, are 



(8) 
(9) 



Mu := u r 



1 



Vu ru r ( u 



(Au) r = 



at r = R, 
at r = R. 



Proof. When fl is a ball, the normal vector to the surface at a point x is n = x/R. 
Then the zth component of (D 2 u)n is given by 



d 

i=i 



and can be rewritten as 



Therefore, 
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Then the projection Pqq takes the tangential component of the above gradient 
vector, and so 

Pan [(D 2 u)n] = graded {R) (Du ■ -| - . 

We know divgogradgQ — Aon by definition. For the ball of radius R, we have 
^dM d (R) — 51 ^S- The operator As is the spherical Laplacian, consisting of the 
angular part of the Laplacian. It satisfies the identity A = + + p-Ag. 

Thus 

div a o-Pao [{D 2 u)n] = A md{R) [Du ■ ~ - ~J 

by noting that Du ■ n = du/dr. The boundary conditions of Proposition [6] then 
simplify to ^ and ([9|, as desired. □ 

The one-dimensional case. The one-dimensional analog of the free plate is the 
free rod, represented by an interval / = [a, b] on the real line. We include its 
boundary conditions for the sake of completeness. 

We may derive the natural boundary conditions from the weak eigenvalue equa- 
tion as before. We obtain as boundary conditions 

b 

= 

a 
b 

= 

a 

and the eigenvalue equation 

Note that these are in fact the one-dimensional analogues of the boundary con- 
ditions and eigenvalue equation obtained in Proposition [6j The computations are 
straightforward integrations by parts and thus omitted. 

The fundamental tone of the free plate in dimensions d > 2 had simple angular 
dependence; the fundamental tone of the free rod under tension can be proved to 
be an odd function. See [71 Chapter 7]. 

Note that we do not have an isoperimctric inequality for the free rod, because 
all connected domains of the same area are now intervals of the same length, and 
are identical up to translation. 



6. The fundamental tone as a function of tension 

Fix the smoothly bounded domain f2. We will estimate how the fundamental 
tone loi — wi(t) depends on the tension parameter r, establishing bounds used in 
the proof of Theorem [T] We will also examine the behavior of uj\ in the extreme 
case as r — > 00. 

First we note that the Rayleigh quotient ([T]) is linear and increasing as a function 
of r. Our eigenvalue ui% (r) is the infimum of Q [u] over u € H 2 (il) with J n udx = 0, 
and thus wi(r) is itself a concave, increasing function of r. 

Next, we will prove wi(r)/r is bounded above and below for all t > 0. Recall 
fix is the fundamental tone of the free membrane. 



12 L. M. CHASMAN 




Lemma 8. For all r > we have 



(10) r/il < Wi(t) < T 



where x = J^xdx/\Q\ is the center of mass ofQ. In particular, when SI is the unit 
ball, 

(11) 77x1 < cji(t) < r(<2 + 2). 



Furthermore, the upper bounds in (10) and ( |ll| hold for all r G K. 
These bounds are illustrated in Figure [2] 

Proof. To establish the upper bound, take the coordinate functions as trial func- 
tions: Uk — Xk — Xk, for k = 1, . . . ,d. Note J n Uk dx — by definition of center of 
mass, so the Uk are valid trial functions. All second derivatives of the Uk are zero, 
so we have 



wi(r) < Q[ufe] 



J n r|L»u fc | 2 da: _ J Q 1 da; 



In u k dx J n (x k - Xk) 2 dx' 

Clearing the denominator and summing over all indices k, we obtain 

oji(t) / |x — x\ 2 dx < r|Sl|d, 



which is the desired upper bound. When SI is the unit ball, note \x\ 2 dx — 
\n\d/(d+2). 

Now we treat the lower bound. Let u G if 2 (SI) with f Q udx — 0. Then 

r , t L \Du\ 2 dx 
Q[u] > J "' ' > r/X! 

by the variational characterization of [i\ . Taking the infimum over all trial functions 
u for the plate yields oj\{t) > t/xi. □ 
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Note that Payne [22! proved linear bounds for eigenvalues of the clamped plate 
under tension. Kawohl, Levine, and Velte [TT] investigated the sums of the first 
d eigenvalues as functions of parameters for the clamped plate under tension and 
compression. 

We can also prove another linear upper bound on u)\, which is just a constant 
plus the lower bound in Lemma [HJ 

Lemma 9. For all t£M, 

UJl < C(fi) + T/il, 

where the value 

L\D 2 v\ 2 dx 

c (n = J " ' 

Jo v dx 

is given explicitly in terms of the fundamental mode v of the free membrane on Q. 

Proof. Let v be a fundamental mode of the membrane with Av = —fJ>iV and 
J^vdx — 0; the membrane boundary condition is du/dn — on dil. Then by 
the variational characterization of eigenvalues, 

wi(r) < Q[v] = C(n) + tQ m [v] = C(O) + r Ml , 

as desired. □ 



6Jj (t) 




r 



Figure 3. The fundamental tone of the disk (solid curve) together 
with the upper bound of Lemma [9] (top dashed line) and the lower 
bound of Lemma [8] (bottom dashed line). 



Infinite tension limit. A plate behaves like a membrane as the flexural rigidity 
tends to zero, that is, as r = (tension/flexural rigidity) tends to infinity. For the 
fundamental tone, that means: 

Corollary 10. For the fundamental tone of the free plate, 

Wi(t) 

> Hi as t — y c*o . 
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Proof. By Lemmas [S] and |9j we have 

wi(t) C 
t r 

Let t — ^ oo . □ 

The eigenfunctions should converge as r — > oo to the eigenfunctions of the free 
membrane problem. Proving this for all eigenfunctions seems to require a singular 
perturbation approach, which has been carried out for the clamped plate in [8], 
but we will not need any such facts for our work. For the convergence of the 
fundamental tone of the clamped plate to the first fixed membrane eigenvalue, see 

HQ. 

Remark. In the limit as r — > 0, we find a relationship between uj\(t)/t and the 
scalar moment of inertia of the region f2; see [7]. 

7. Summary of Bessel Function facts 

The radial part of the fundamental tone of the unit ball is a linear combination 
of Bessel functions. Before we begin constructing our trial functions, we need to 
gather some results established in the companion paper [BJ. 

The ultraspherical Bessel functions ji (z) of the first kind are defined in terms of 
the Bessel functions of the first kind, J u {z), as follows: 

ji(z) = z~ s J s+ i(z) 
d-2 

with s = . 

2 

and solve the ultraspherical Bessel equation 

(12) z 2 w" + (d - l)zw' + (z 2 -l(l + d-2))w = 0.. 

Note that this notation suppresses the dependence of the ji functions on the di- 
mension d; we assume dimension d > 2 is fixed. Ultraspherical modified Bessel 
functions ii(z) of the first kind are defined analogously, with 

it(z) = z~ s I s+ i(z) 

solving the modified ultraspherical Bessel equation 



(13) z V + (d - l)zw' - [z 2 + 1(1 + d - 2)) w 

L 
(14) 



= 0.. 

Ultraspherical Bessel functions satisfy the following recurrence relations: 
d - 2 + 21 





= 3l-i(z) 


+ .ii+i{z) 




= l -m- 

z 


- 3i+i{z) 




= k-i(z) 


-ii+i{z) 


i'i(z) 


= -h{z) - 


1- k+i{z) 



z 

(15) 

(16) d ^ 2 + 2! 

z 

(17) 

z 

Ultraspherical Bessel functions and their derivatives may be expressed by con- 
verging series. The first few terms in these expansions may be used to bound the 
Bessel functions and their derivatives; we will need such bounds for the second 
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derivatives of fx and i\. Let dk denote the coefficients of the series expansion for 
i"(z), so that 



Ji (*) = X^" 1 )*** 2 *" 1 and ii'(z) = J2 dkz 2 "- 1 

k=l k=l 

by the series expansions of ji and i\ in [SI P- 5] , where 

i _ 2fc + 1 2 1 -2 fc - d /2 

k (k - l)W(k + 1 + d/2) 

Lemma 11. [HI Lemma 10] We /lave i/ie following bounds: 

-d x z + d 2 z 3 > jj'(z) for all z £ [o, + 2)/(d + 5) 



6 

d\z H — c?2 2;3 > *i ( z ) / or oZZ 2 S 

5 



0,^3 



While proofs are provided in the companion paper, these bounds and those listed 
below in Lemma [12] can also all be demonstrated numerically in Mathematica or 
Maple for any given dimension. 

We will also be using additional facts about the signs of certain Bessel functions 
and derivatives. These were proven in [6] and are collected below. We write doo for 
the first nontrivial zero of j[. 

Lemma 12. 6j Lemmas 5 through 9] We have the following: 

(1) For I = 1, ... ,5, we have ji > on (0, cioo]. 

(2) We have j[ > on (0, a^). 

(3) We have j' 2 > on (0, aoo]. 

(4) We have j'{ < on (0, o<x,]. 

(5) We have jf^ > on (0, doo]. 

We are now ready to begin proving Theorem [TJ 

8. Trial functions 

In this and the next sections, we establish the lemmas which allow us to prove 
Theorem [TJ Among all regions fl of a fixed volume, when r > the fundamental 
tone of the free plate is maximal for a ball. That is, 

(18) Wi(f2) < u!i(fl*), with equality if and only if SI is a ball. 

For simplicity, as we prepare to prove Theorem [TJ we will write ui instead of u>i 
for the fundamental tone of the free plate with shape Q; the fundamental tone of the 
unit ball will be denoted by u>* . When dependence on the region Q and the tension 
r need be made explicit, we write ui(t, fl) for the fundamental tone and Q T .n for 
the Rayleigh quotient. The tension parameter r > throughout the remainder of 
this paper. 

We begin with the assumption that our domain SI has the same volume as the 
unit ball M d ; this is justified by the scaling argument in Lemma 24 

In [51 Theorem 2] , we identified the fundamental mode of the unit ball, written 
in spherical coordinates as: 

(19) m = R{r)Y x := (ji(ar) + ^h(br)^Y u 
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where Y% is any of the d spherical harmonics of order 1, ji(z) and i\{z) are ultras- 
pherical Bessel functions, and a and b are positive constants satisfying the conditions 
a 2 b 2 = to and b 2 — a 2 = r. Note that we may take our spherical harmonics Y\ to 
be the Xi/r for i = 1 . . . , d, where Xi the ith coordinate function. 

Finally, recall that we took to be the first nontrivial zero of j[(z); by the 
proof of Theorem 3] , we have a < . Note that the fundamental tone of the 
free membrane with shape M d is given by p,\ = a 2 ^ . 

We are now able to choose our trial functions. Inspired by Weinberger's proof 
for the membrane [3 2) , we choose appropriate trial functions from the fundamental 
modes of the unit ball. In the following lemmas, we take 

R( r ) = 3i( ar ) +jii(br) 

to be the radial part of the fundamental mode of the unit ball. Recall a and b are 
positive constants determined by r and the boundary conditions, as in the proof 
of Theorem 3 in [BJ. The constant 7 is positive and determined by a, r, and the 
boundary conditions as folllows: 

(20) 7 •= > 0. 

Recall also that R(r) > on (0, 1] and R'{1) > 0. 

Lemma 13. (Trial functions) Let the radial function p be given by the function R, 
extended linearly. That is, 



p{r) 



R(r) when < r < 1, 

R(l) + (r-l)R'(l) whenr>l. 



After translating Q suitably, the functions = Xkp{r)/r, for k = l,...,d, are 
valid trial functions for the fundamental tone. 

Proof. To be valid trial functions, the Uk must be in H 2 {yi). Because is bounded 
in R d , the only possible issue would be a singularity at the origin. The series 
expansions given in [SJ p. 5] for ji and i\ give us that R(r)jr approaches a con- 
stant as r — > 0. Thus, Uk € H 2 (tt) as desired. The trial functions must also be 
perpendicular to the constant function, and so we will need 

E^th.dx = for k=l,...,d. 
r 

We use the Brouwer Fixed Point Theorem to translate our region so that the 
above conditions are guaranteed; here again we follow Weinberger [32]. Write 
x = (x\, . . . , Xd) and consider the vector field 

X(v)= ( p{ } X ~ V }\ x-v)dx. 



The vector field X is continuous by construction. For any vector v along the 
boundary of the convex hull of f2, the vector field X(y) is inward-pointing, because 
p > and the entire region ft lies in a half-space to one side of v. Thus by the 
Brouwer Fixed Point Theorem, our vector field X vanishes at some v in the convex 
hull of f2. If we first translate fi by v, then we have X(0) = J Q p(r)x/rdx = 0. 
This gives us J n Uk dx — 0, as desired. □ 

We will need one further fact about our radial function p. 
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Lemma 14. (Concavity) The function p"(r) < for r 6 [0, 1], with equality only 
at the endpoints. 

Proof. First note that on [0, 1], the function p = R. We see 

R"(r)=a 2 rf(ar) +1 bX(br), 

which is zero at r = because the individual Bessel derivatives vanish there, by 
the series expansions for the Bessel ji(z) and ii{z) in [6]. At r = 1, the function R" 
vanishes because of the boundary condition Mu = 0. 
The fourth derivative of R is given by 

R""{r) = a i j'{"{ar)+ 1 bH'{"{br). 

Because all derivatives of ii(z) are positive when z > 0, the second term above 
is positive on (0, oo). Lemma 9 in [6] states that that j""(z) is positive on (0, a<x,]. 
Thus R""{r) > on (0, 1], and so R"{r) is a strictly convex function on [0, 1]. Since 
R" = at r = and r = 1, the function R" must be negative on the interior of the 
interval [0,1]. □ 

We now bound our fundamental tone above by a quotient of integrals whose 
integrands are radial functions. The numerator will be quite complicated, so we 
write 

; = (P ) + — — [p-rp) +t{p) + — - 2 — p. 
We will also need the following calculus facts: 
Fact 15. [7, Appendix] We have the sums 

EKI 2 = p 2 



fc=l 

(/ 



J\2 



k=l 

Y j \D^ = { P "f + ^^ 1 {p-rp'f 



k=l 



We may now use the trial functions to bound our fundamental tone by a quotient 
of integra;s. 



Lemma 16. (Using the trial functions) For any VL, translated as in Lemma 13, we 
have 



(21) w < 



In N IpI dx 



In P 2 dx 
with equality if Q = Q* . 

Proof. For dehned as in Lemma [13] we have 

, J n \D 2 u k \ 2 + r\Du k \ 2 dx 
U ~ Q[Uk] = f ImP da: ' 
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from the Rayleigh-Ritz characterization. We have equality when fl = f2* because 
the Uk are the eigenfunctions for the ball associated with the fundamental tone; see 
[B]. Multiplying both sides by J n \uk\ 2 dx and summing over all k, we obtain 

d „ d d 

(22) w / J2\u k \ 2 dx< / V|^ 2 Mfe | 2 +rV|^ fc | 2 dx 

Jn ,. , in,. , 



fe=i 



again with equality if = fT. 

By these sums in Fact [l5j we see inequality ( 22 1 becomes 

P 2 dx< ^(W +' i ^^{ P -rp'f + r(p'f + ^ 

once more with equality if ft is the ball fl* . Dividing both sides by J s ~ l p 2 dx, we 
obtain " □ 

9. Partial monotonicity of the integrands 



We want to show the quotient (21 ) in Lemma 16 has a sort of monotonicity with 
respect to the region f2, and so we examine the integrands of the numerator and 
denominator separately. The case of the denominator is much simpler; the partial 
monotonicity of the integrand of the numerator is much more difficult, and requires 
several lemmas. 

We begin with the denominator. 

Lemma 17. (Monotonicity in the denominator) The function p{r) 2 is strictly in- 
creasing. 



Proof. Differentiating, we see 

j[(ar) + yi^ibr) when < r < 1, 
R'(l) when r > 1. 



P'(r) 



Obviously i'i(br) > 0. Because we have a < from the proof of Theorem 3 in 
[5J, the function j[(ar) is positive on [0, 1]. Thus p'(r) is positive everywhere, and 
p (and therefore p 2 ) is an increasing function. □ 

We do not need to prove the integrand of the numerator is strictly decreasing; a 
weaker "partial monotonicity" condition is sufficient. We will say a function F is 
partially monotonic for Q if it satisfies 

(23) F(x) > F(y) for all x € n and y O. 

Lemma 18. (Partial monotonicity in the numerator) The function 

N[P] = {P"f + 3 ~^ }1 (P rp>) 2 + r Up') 2 + (d - 1)£ 



satisfies condition ( 23 ) for Q the unit ball. 



Proof. Given that p" < on (0, 1) and equals zero elsewhere by Lemma 14 the 
function (p") 2 satisfies condition (23 1 for the unit ball. The derivative of the func- 
tion t(p') 2 with respect to r is 2rp' p" , and hence negative on (0, r) and zero every- 
where else. Thus t{p') 2 is a decreasing function of r. It remains to show that the 
remaining term 
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is also a decreasing function of r. Differentiating, we see 

h '( r ) = ^P~ r P') C^(p-rp')+3p" +rp 
Now, p — rp' = at r = and 



dr 



(p-rp') = -rp", 



so by Lemma f 4 (p — rp') is positive on (0, oo) and vanishes at zero. Thus in order 
for h{r) to be decreasing, we must have 

6 



(24) 



(p - rp) + 'dp" + rp > 0. 



Let A r p := p" — (d — l)r (p — rp'). Recall from the Bessel equations (12) 
and Jl3l) that 



(25) 



A r ji(ar) = —a 2 ji(ar) and A r ii(br) — b 2 ii(br). 



Then on the interval [0,1], 
4(> - rp') + 3p" + rp = ~{p - rp') + 3 ( A r p + - 



(p - rp') + rp 



= 3 ^ r t ~ (P - rp) + 3( - a 2 ji(ar) + jb 2 ii(brfj + rp, 

with the last equality by (25 ). Considering the first term of the last line above, we 
see by (15) and (17), 

-\{P~ rp') = ^(arj 2 (ar) - brji 2 (br)^j 

= (^(jiM+JsM) +lb 2 (i 3 (br)-i 1 (br)^ 

with the second equality by (14), (16), and simplifying. 

Therefore our quantity of interest in ( 24 ) can be bounded below in terms of ji 's 
and z;'s: 

^{p-rp')+3p" + t P 



3a 2 
d + 2 



+ 7 t + 



3b 



d + 2 



ji(ar) + 
i\(br 



3a 2 (d+l) . 



d+2 
36 2 (d+l) . 



■7" 



d + 2 



J3(ar) 

i 3 (br). 



3a 2 (d+l) . . , 36 2 (d+l). , 



d + 2 



3a 2 
d + 2 



+ 7 r + 



d + 2 
3b 2 
d + 2 



Ji(ar) 



with the inequality by ji(ar) < ii(ar) < i;(6r), since r > and so a < 6. 

The function «3 is everywhere nonnegative and the constant 7 is positive. We 
have a < aoo, so ar < 0,00 on [0, 1] and hence the functions j3(ar) and ji(ar) are 
positive on [0, 1] by Lemma 1 of [BJ. The remaining factor is positive for all r > by 
Lemmas 19 and[20|(to follow), thus establishing (24) and completing the proof. □ 
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We establish the positivity of the remaining factor first for those r values such 
that t > 9/(d + 5); the proof for smaller r values is more complicated and is treated 
in another lemma. 



Lemma 19. (Large t) We have 
for allr> 9/(d+5). 

Proof. We use the bounds we established for oj(t) in Section [4] 

Recall that the first free membrane eigenvalue for the ball is fj,* = a 2 ^. Lemma [8] 
and Proposition of Lorch and Szego (Q2], but see the statement in Prop 4]) 
together give (d + 2)r > to* > rd. Because uo* 
inequalities relating r and a: 

„4 



a 4 + a 2 r [SI Prop 2], we obtain 



a — a z a 4 

with the upper bound holding only if a 2 < d. 
Using the lower bound, we see 



3a 2 



> 



a 



2 -a 2 



3a 2 



d + 2 d+2 



a 4 (d + 5) 



d+2 
3a 2 (d + 2) 



(d + 2)(d + 2-a 2 ) 



which is nonnegative whenever a 2 > 3(d + 2)/(d + 5). When a 2 < 3(d + 2)/(d+5), 
we have 

3a 2 9 



d + 2 



> T - 



by our choice of r. 

Lemma 20. (Small t) We have 

/or a;/0<T<9/(d + 5). 



d + 2 



+ 7 



d + 5 



36 2 
d + 2 



>0, 



□ 



> 



Proof. The proof will proceed as follows. For < r < 9/(d+ 5), we restate the 
desired inequality (28) as a condition on 7, (30 1. We then use properties of Bessel 
functions to establish a lower bound on 7 in terms of a rational function of a; we 
then show this function satisfies ( 30 ) . We will need to treat the cases of d > 3 and 
d = 2 separately, because the two-dimensional case requires better bounds than we 
can derive for general d. 
First note that b 2 = 



(29) 



r, so the inequality (28) is equivalent to 
3a 2 (1 



T > 



7) 

(d+ 5)7+ (d + 2)' 



Using the lower bound on r in (|27|), we see that the above will hold if 
(30) 7 > 



3(d + 2)-a 2 (d + 5) 



(3 + a 2 )(d + 2) 



=: 7 
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We need only show that ( 30 ) holds for all < r < 9/(d + 5). We will use Taylor 
polynomial estimates to bound 7 below by a rational function. From Lemma 6 of 
[BJ, we have 



i"{z) < —diz + d 2 z 3 
6 

i"{z) < diz + -d 2 z 3 
5 



[(W3(d + 2)/(d + 5)], 
on [0,V3],. 



These bounds apply to z = ar and z = br respectively, when r € [0, 1], as we show 
below by obtaining bounds on a 2 and b 2 . 

To derive our bound on a 2 , we note that the lower bound of (27 1 together with 
our assumption r < 9/(rf + 5) implies 



d + 2 - a 1 



< 



d + 5' 



so that 



(d + 5)a 4 + 9a 2 -9(d+2) < 0. 



The left-hand side is increasing with respect to a 2 and equals zero when a 2 = 
3(d + 2)/(d + 5). Hence a 2 < 3(d + 2)/(d + 5) and the bound on j'({z) holds for 
z = ar when r < 9/(d + 5). We use these to obtain a further bound: 



> r - 



d+5 
,4 



> 



d + 2 - a 2 



d + 2 - a 2 



9 

2 3(rf + 2) 
2 d + 5 



and so we have a 2 < d. 

To bound 6 2 we use b 2 = ■ 



t and obtain 



6 2 = 

and so b 2 < 3. 

We also have, from (27 1, 

(31) 



r < 



3(d + 2) 9 
d + 5 + d + 5 



= 3, 



da 2 , 2 (d + 2)a 2 
^ „ > 6 > 



d + 2 - a 2 



with the upper bound holding in this regime because a 2 < d. 
We also need the following binomial estimate: 



(32) 



1 - -x < (1 - xf' 2 
2 v ' 



for < x < 1. 
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Using these bounds, we see 

= -« 2 jT(«) 
7 b 2 i'{(b) 

a 2 (dia — d 2 a 3 ) 



> 



> 



b 2 (d 1 b+{6/5)d 2 b 3 ) 



by definition ( 20 ) 



by Lemma 6 of [6] 



a 3 (d\ — d 2 a 2 ) 



da 2 
d—a 2 



3/2 



(di 

3/2 



> 1- 



+ (6/5)rf 2 ^) 

(rf-a 2 )(l- Cl a 2 ) 
(d-a 2 + (6/5)cida 2 
Sa 2 ^ (d-a 2 )(6(d + 4)~5a 2 ) 
~2d 



by (31) 



writing ci = d 2 jd\ 



6(d + 4) 



by (32), 



(6d(d + 4) - 24a 2 ) 

noting that a 2 /d < 1 and a 2 < 3(d + 2) /(d + 5). 
Thus we have 7 — 7* > if 

3a 2 \ (d - a 2 )(6(d + 4) - 5a 2 ) 3(d + 2) - a 2 (d + 5) 



1 



2d 



>0, 



6d(d + 4)-24a 2 (3 + a 2 )(d + 2) 

or, clearing the denominators and writing x = a 2 , if 

(2d-3x)(d-x)(6(d+4)-5x^(3+x)(d+2)-2d(6d(d+4)-24x^ (3(d+2)-x(d+5)) > 0. 

The above polynomial is fourth degree in each of d and x and has the root x = 0; 
because we are only interested in its behavior for x 6 (0, 3(d + 2)/(d + 5)), we may 
divide by x and work to show the resulting polynomial 



P(x, d) = 24d 4 + 60d 3 - 120d 2 - 432d - 40d 3 x 
+ 43d 2 x 2 + 113dx 2 + 54a; 2 - 15dx 3 



119d 2 x 
-30x 3 



6dx + 432a; 



is nonnegative for x G (0, 3(d + 2)/(d + 5). This claim is addressed in Lemma 21 
for d > 3. 

For d = 2, the function P(x, 2) is negative on most of our interval of interest 
[0,12/7], and so we must improve our lower bound on 7. The derivation follows 



that of inequality ([27]) in the proof of Lemma 19 as follows 
By Lemma |8j u>* > aLr, where 



1.84118 is the first zero of Jx{z). By 



Proposition 2 of [5] we have lu* = a 4 



a r, giving us 
a 4 



Using b 2 



T < 

r, we obtain also a bound on b 2 : 



b 2 < 



2 2 



Proceeding as before, we deduce 



7> 1 



3a 2 
2al 



(ag -a 2 )(36-5a 2 ) 
360^ + (60^ - 36)a 



with the last again from ( 32 ) . So 7 — 7* > if 
3a 2 \ (a 2 -a 2 V36- 



5a 2 ) 



12 - 7a 2 
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or, setting x = a , if the fourth degree polynomial 

Q(x) =(l- (al - x)(36 - 5x)(12 + Ax) - (zGal + (6a 2 , - 36)z) (12 - 7x) 

is positive on [0,12/7]. This positivity follows from Lemma 22 completing our 
proof. □ 

The next two lemmas regarding the polynomials P and Q allow us to complete 
the proof of Lemma [20| 

Lemma 21. The polynomial 

P(x, d) = 24d 4 + 60d 3 - 120d 2 - 432d - 40d 3 x - 119d 2 x - 6dx + 432a; 

+ 43d 2 x 2 + 113dx 2 + 54x 2 - I5dx 3 - 30x 3 

is nonnegative for all x € (0, 3(d + 2)/(d + 5)) and integers d > 3. 

Proof. First note that 3(d + 2)/(d + 5) < 3. We bound P below on the interval 
x £ [0, 3] by taking x = 3 in terms with negative coefficients and taking x = in 
terms with positive coefficients, obtaining 

P{x, d) > 24d 4 - 60d 3 - A77d 2 - 855d - 810 =: g(d). 

The highest order term is 24d 4 , and so g is ultimately positive and increasing in d. 
Note also that 

g'{d) = 96d 3 - 180d 2 - 954d - 855 

g"(d) = 288d 2 - 360rf- 954. 

The function g 1 ' '(d) is a quadratic polynomial with positive leading coefficient and 
roots at d « —1.30 and 2.55; thus g'(d) is increasing for all d > 3. We see that 
g'{5) = 1875, so g is increasing for all d > 5. Finally, g(7) = 6876, so for all d > 7 
we have g(d) > and hence P(x, d) > for all d > 7 and x <E [0, 3]. 

For d = 3,4,5,6, we look at the polynomials Pd(x) = P(x,d) directly to show 
that Pd(x) > on [0, 3(d + 2)/(d + 5)]. Each Pd is a cubic polynomial in x; its first 
derivative P d (x) is quadratic and so the critical points of Pd(x) can all be found 
exactly. 

For d = 4, 5, and 6, direct calculations show P' d < on [0, 3(d + 2)/(d + 5)] and 
P d {3(d + 2)/(d + 5)) > 0, so P d {x) > on [0, 3(d + 2)/(d + 5)]. 

For d = 3, our interval of interest is [0, 15/8]. We have a critical point 
c fa 1.4 e [0, 15/8], with P^ < on [0,c] and P3 > on [c, 15/8]. The critical value 
P 3 (c) « 79 is positive, so P%(x) > on the desired interval [0, 15/8]. □ 

Lemma 22. The polynomial 

Q(x) = (l- ~) (al - x)(36 - &r)(12 + 4x) - (36a 2 , + (6a 2 , - 36)x) (12 - 7x) 
is positive on [0, 12/7]. 

Proof. As in previous cases, x = is a root of this polynomial, so we examine 
(?(x) := Q(x)/x. The derivative g'(x) is a quadratic polynomial, so its roots can be 
found exactly. We see that g has a critical point c pa 1.4 in [0, 12/7], with g' < on 
[0, c] and 3' > on [c, 12/7]. The critical value g(c) sa 177.8 is positive, so g > on 
[0,12/7]. □ 
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10. Proof of the isoperimetric inequality 

Now that we have established the desired monotonicity of our quotient, we need 
two more lemmas before we can prove the isoperimetric inequality for the free 
plate under tension. Our next lemma is a simple observation about integrals of 
monotone and partially monotone functions, which is a special case of more general 
rearrangement inequalities (see [El Chapter 3]). 

Lemma 23. For any radial function function F(r) that satisfies the partial mono- 



tonicity condition (231 for ft* , 



Fdx < Fdx 
n Jn* 



with equality if and only iffl — fl*. For any strictly increasing radial function F(r), 



Fdx> Fdx 
n Jn* 



with equality if and only Q = Q* 



Proof. Note that = with = Suppose F satisfies (|23) for ft* 

The result follows from decomposing the domain: 



Fdx= Fdx+ Fdx 

n Jnnn* Jn\n* 



< Fdx+ sup |F(x)||fi\fi*| 
Jnnn* xen\n* 

< Fdx+ inf |F(a;)||n*\n| since F satisfies ([23]). 



< / Fdx+ Fdx 
Fdx. 

n* 

Note that if > 0, either the second inequality or the third is strict by the 



strict inequality in ( 23 1 . If F is increasing, then apply the first part of the Lemma 



to the function —F. □ 



The final lemma describes how the eigenvalues change with the dilation of the 
region, and is used in the proof of the theorem to show we need only consider f2 
with volume equal to that of the unit ball. We will use the notation sil := {x € 
R d : x/s G il} for s > 0. 

Lemma 24. (Scaling) For all s > 0, we have 



ui(t,CI) = s 4 U](s- 2 t,sQ). 
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Proof. For any u £ H 2 {Vl) with J^udx = 0, let u(x) — u(x/s). Then u is a valid 
trial function on s^l and so 

j sn \p 2 u\ 2 + s - 2 t\du\ 2 dx 

_ IsQ \ s ~ 2 ( d2u )( x / s )\ 2 + s- 2 T\s- l {Du){x/s)\ 2 dx 
J sn u(x/s) 2 dx 



s- 4+d J n \D 2 u\ 2 +r\Du\ 2 dy 
s d J n u 2 dy 

s~ 4 <9r,nM- 



taking y = x/s, 



Now the lemma follows from the variational characterization of the fundamental 
tone. □ 

We can now prove our main result. 



Proof of Theorem^ Once we have established inequality (18) for all regions il of 
volume equal to that of the unit ball and all r > 0, we obtain ( 18 ) for regions of 
arbitrary volume, since 

lj(t, n) = s i Lo{ s - 2 T 1 so) < s 4 u}( s - 2 t, s n*) = lj(t, n*), 

for all s > by Lemma [24] 

Thus it suffices to prove the theorem for Q with volume equal to that of the unit 
ball, so that 0* is the unit ball. We may also translate 51 as in Lemma [13J which 
leaves the fundamental tone unchanged. Then, 

L N\p\ dx , 

u < r , , by Lemma HE 
Jn P dx 

< in* N \p]dx Lemmas fTTl fl8l and [23l 

J n ,p 2 dx i—ii—i i—i 

= w , 



by applying the equality condition in Lemma 16 Finally, if equality holds, then f2 
must be a ball, by the equality statement in Lemma [23) □ 



11. Further Directions 

The isperimetric problem for the free plate considered in this paper can be gen- 
eralized in several different directions: considering the case where the material 
property Poisson's ratio is nonzero, investigating a stronger inequality involving 
the harmonic mean of eigenvalues, and considering the problem on curved spaces. 

Poisson's Ratio. One generalization of the free plate problem is to account for 
Poisson's Ratio, a property of the material of the plate that describes how a rec- 
tangle of the material stretches or shrinks in one direction when stretched along 
the perpendicular direction. Our Rayleigh quotient and work so far all hold for a 
material where Poisson's Ratio is zero. Most real- world materials have a E [0, 1/2], 
although there exist some materials with negative Poisson's Ratio. 
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We will assume a G [0, 1) in order to be assured of coercivity of the generalized 
Rayleigh quotient, given by 

_ f n (l-a)\D^ + a(Aur + r\Du\"dx 

(33) " JoH^ • 

This quotient reduces to our previous quotient ([I]) when a = 0. Following our 
earlier derivation, we obtain the same eigenvalue equation 

AAu — tAu = LOU, 

along with new natural boundary conditions on dfl, which reduce to the old ones 
when a = 0. 



(l-a)—+aAu 
an 1 

- (1 - a)div dQ (P 9n [{D 2 u)n}) - ^ 
On y 1 11 On 



= 0, 

an 

= 0. 

an 



The generalization to nonzero o does not change the eigenvalue equation and 
hence the general form of solutions is preserved. However, the change in the 
Rayleigh quotient affects the proof of Theorem 3 of [S] , which identified the funda- 
mental mode of the ball. This in turn affects the proof of the isoperimetric inequality 
in this paper. We can no longer complete the square in the Rayleigh quotient as in 
Theorem 3] to show the fundamental mode of the ball corresponds to I = or 1, 
although for some values of a we can adapt the proof to show the lowest eigenvalue 
corresponding to I = 1 is lower than that for / = 0. 



Harmonic mean of low eigenvalues. In two dimensions, Szego was able to prove 
a stronger statement of the Szego- Weinberger inequality using conformal mappings 
Pr?l 128] . Specifically, he proved that the sum of reciprocals 

1 1 

h — 

Mi A*2 

is minimal for a disk. In other words, the harmonic mean of and /i 2 is maximal 
for the disk. Our investigation in Chapter 3] with the moment of inertia suggests 
a similar result for the free plate, since the moment of inertia is minimal for a ball. 
That is, for the free plate, we conjecture 

- V — > 1 



Curved spaces. We have taken our region f2 to be in Euclidean space M. d , but 
we could consider the same eigenvalue problem on a region in spaces of constant 
curvature: the sphere and hyperbolic space. Other eigenvalue inequalities have 
been proven in these spaces [jQ. In particular, the Szego- Weinberger inequality was 
proved for domains on the sphere by Ashbaugh and Benguria [2 . Another direction 
of generalization would be Hersch-type bounds for metrics on the whole sphere or 
torus; see [13] . 
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